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Abstract. We report upon the present status of global fits to Cabibbo-Kobayashi-Maskawa matrix. 
PACS. XX.XX.XX No PACS code given 



1 Introduction 

The three-family Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix is a key element of the Standard 
Model (SM). The nine complex CKM elements are com- 
pletely specified by three mixing angles and one phase 
that is responsible for CP violation in the SM. Measuring 
the CKM matrix elements in various ways provides consis- 
tency tests of the matrix elements itself and with unitarity. 
Any significant inconsistency with the SM would indicate 
the presence of new physics. 

A convenient parameterization of the CKM matrix is 
the Wolfenstein approximation 1 , which to order 0(A 3 ) 
is given by: 
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where A = 0.2241 ± 0.0033 is the best-known parameter 
measured in semileptonic K decays, A = 0.82 is deter- 
mined from semileptonic B decays to charmed particles 
with an accuracy of c± 6% and p = p • (1 — A 2 /2) and 
fj = rj - (1 — A 2 /2) are least-known. 

The unitarity of the CKM matrix yields six trian- 
gular relations of which V ud V* b + V cd V* b + V td V* h = 
is well-suited for experimental tests. In order to deter- 
mine the apex of the unitarity triangle (p, fj) presently 
eight measurements are used as input, the B semileptonic 
branching fractions B(B — > X c tv\ B(B — > X u tv\ and 
B(B — > piv\ the normalized B — > D*£v rate at zero re- 
coil, T(l)\V u b\ 2 , the B® and B® oscillations frequencies 
AmB d and Z\ra# s , the parameter \^k\ that specifies CP 
violation in the K°K° system, as well as sin 2/3 which is 
measured in CP asymmetries of charmonium Kg (K L ) 
final states. Though many of these measurements them- 
selves are rather precise their translation to the p — fj plane 
is affected by large non-gaussian theoretical uncertainties. 
Various approaches, which treat theoretical errors in dif- 
ferent ways, can be found in the literature [2,3,4,5,6]. 



2 The Scan Method 

The scan method is an unbiased procedure for extract- 
ing A,p,fj from measurements. We select observables that 
allow us to factorize their predictions in terms of theo- 
retical quantities Ti that have an a priori unknown (and 
likely non-gaussian) error distribution (Ai), other observ- 
ables, and the CKM dependence expressed as functions 
of Wolfenstein parameters. As an example, consider the 
charmless semileptonic branching fraction for B — > piv, 
which is predicted to be B(B — > ptv) = |^6| 2 • P p t v • 
tb, where r B o is the B° lifetime and T p t v is the reduced 
rate affected by non-gaussian uncertainties. This analysis 
treats eleven theoretical parameters with non-gaussian er- 
rors, the reduced inclusive semileptonic rates Px u tv and 
Px c tv, the form factor for B — > D*tv at zero recoil, Tb* (1)? 
the bag factors of the K° and B° systems, Bk andl^g, 
the B° decay constant f B , £ 2 = f B J 'f Bd B B J 'B Bd and 
the QCD parameters 7/1,772,773 and 77^. 

We perform a x 2 minimization based on a frequentist 
approach by selecting a specific value for each Ti within 
the allowed range (called a model). We perform individ- 
ual fits for many models scanning over the allowed non- 
gaussian ranges of the Ti parameter space. The QCD pa- 
rameters are not scanned; their small errors are treated in 
the x 2 as gaussian. For theoretical quantities calculated on 
the lattice, which have gaussian errors (Bk, Bb, /b and £) 
we add specific x 2 terms. To account for correlations be- 
tween observables that occur in more than one prediction, 
such as the masses of the t-quark, c-quark, and VF-boson, 
B hadron lifetimes, B hadron production fractions and A, 
we include additional terms in the y 2 function. 

We consider a model to be consistent with the data if 
the fit probability yields P(x 2 ) > 5%. We determine the 
best estimate for each of the 17 fit parameters and plot a 
95% confidence level (C.L.) contour in the p — fj plane. We 
overlay the p — fj contours of all accepted fits. In order to 
study correlations among the Ti and constraints the data 
impose we perform global fits with non-gaussian theory 
errors scanned over a ±5 A wide range (see section 5). 
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Table 1. Measurement inputs used in x 2 minimization 



Observable 



Value 



Comment 



B(B X u iv) (2.03 ± 0.22 exp ± 0.31 th ) x 10" 

B(B -> X u tv) (1.71 ± 0.48 e * p ± 0.21^) x 10" 

B(B X c iv) 0.1070 ± 0.0028 

B{B -> Xdv) 0.1042 ± 0.0026 

B(B -> piv) (2.68 ± 0A3 exp ± 0.5th) x 10" 3 

|y cfe |F(l) 0.0388 ± 0.005 ± 0.009 

Z\m Sd (0.503 ± 0.006)ps _1 

Am Bs > 14.4 p S - 1 @95%C.L. 

\e K \ (2.271 ±0.017) x 10" 3 

sin 2/3 0.731 ±0.055 

A 0.2241 ± 0.0033 
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Fig. 1. Results of the global fit in the p — r\ plane. 



Table 2. Theoretical parameter with non-gaussian errors 



0.87 < T D *{1) < 0.95 
12.0 < r(piu) < 22.2 ps" 1 



38.0 < r(c£is) < 41.5 ps" 1 
54.8 < f{u£v) < 79.6 ps -1 



0.72 <B K < 1.0 (t Bk 
211 < f Bd \[~EhT d < 235 MeV a 
1.18 < g < 1.30 (7$ = 



'b b 
0.04 



0.06 (gaussian) 

33 MeV (gaussian) 
gaussian) 



0.54 < j] B < 0.56 
0.564 < r] 2 < 0.584 



1.0 < 771 < 1.64 
0.43 < 773 < 0.51 




Fig. 2. Fit results in S — p plane from a 



ct>K% 



2.1 Treatment of ArriB s 



Since B®B® oscillations have not been observed yet, a 
lower limit on ArriB s at 95% C.L. has been determined 
by combining analyses of different experiments using the 
amplitude method 8 j. To incorporate ArriB s into the x 2 
function, we use a new approach that is based upon the 
significance of a ArriB s measurement |J2|: 



2w)e 



-±(Am s cr t ) 
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where N is the sample size, Jb s is the B s purity, w is 
the mistag fraction, and a t is the resolution. Substituting 

C for y^/s s (l — 2w) and interpreting S as the number 
of standard deviations by which ArriB s differs from zero, 
S = AvtlbJ o Am Bs •> we ma y define a contribution to the 
2 from the ArriB s measurements as: 
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where A is the best estimate according to experiment. 
The values of (Z\, C 2 , a t ) are chosen to give a minimum at 
17 ps -1 , and a P(x 2 ) = 5% at Am Bs = 14.4 ps -1 . In the 
region of small x 2 ? this function exhibits similar general 
features as that used in our previous global fits ^01 5 while 
it does not suffer from numerical instabilities arising from 
multiple minima. The two functions deviate at large values 
of x 2 , where in any case poor fits result. 



Table 3. Results of 95% C.L. range for p,rj,a and 7 from the 
global fits shown in figure □ For comparison results from RFIT 
and the Bayesian method are also given. 



parameter 


Scan method 


RFIT 9 


Bayesian 9 


P 


-0.13 to 0.40 


0.091 to 0.317 


0.137-0.295 


fj 


0.22 to 0.48 


0.273 to 0.408 


0.295-0.409 


a 


50.4° to 126.6° 






7 


34.4° to 91.7° 


42.1° to 75.7° 


47.0° to 70.0° 



3 Results of the global Fit 

Figure shows the result of scanning all Ti simultane- 
ously within ±1Z\ of their allowed range except for the 
QCD parameters. We have used the input measurements 
summarized in tableland ranges for the T$ listed in ta- 
ble The black points represent the best estimates of 
(p, fj) for each model that is consistent with the data. 
The grey region shows the overlay of all corresponding 
95% C.L. p — fj contours. For reference, the light ellipse 
depicts a typical contour. To guide the eye the 95% C.L. 
bounds on \ € k\, ^ m B d and sin 2j3 as well as the 

lower bound on ArriB s are also plotted. From these fits 
we derive 95% C.L. ranges for p,fj,a and 7 that are listed 
in table 01 For comparison, recent results from two other 
global fits (RFIT 0|, Bayesian fit jHj) are also shown. 

Using the same source of inputs, several differences 
exist between the scan method and the other two ap- 
proaches. First, we scan separately over the inputs of ex- 
clusive and inclusive b — > utv and b — > civ measurements. 
Second, we use a different approach to incorporate AmB s - 
While in the Bayesian method theoretical quantities are 
parameterized in terms of gaussian and uniform distribu- 
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tions, we make no assumptions about their shape. Thus, 
the Bayesian fits tend to produce a smaller region in the 
p — fj plane and are more sensitive to fluctuations than 
corresponding fits in the scan method. In RFIT, the p — fj 
plane is scanned to find a solution in the theoretical pa- 
rameter space. Since in RFIT a central region with equal 
likelihood is determined, it is not possible to give probabil- 
ities for individual points. In contrast, in the scan method 
individual contours have a statistical meaning, with the 
center point yielding the highest probability. Since the 
mapping of the theory parameters to the p — fj plane is not 
one-to-one, it is possible in the scan method to track which 
values of (p, fj) are preferred by the theory parameters. 

4 Search for New Physics 

The decay B — > c/)Kg that proceeds via a b — > sss pen- 
guin loop is expected measure sin 2(3 in the SM to within 
<~4%. New physics contributions, however, may introduce 
a new phase S that may change the CP asymmetry a^o 
significantly from aj^ K o. The BA BA Rj Belle average of 
S(f)K — —0.39 ±0.41 has been updated this summer yield- 
ing S^o = -0.14 ± 0.33 Q3]. The deviation from sin 2(3 
has remained at ~ 2.6<r. In our global fit we introduce 
a new phase S . Figure shows the overlay of all result- 
ing contours in the S — p plane that have acceptable fit 
probabilities. Presently, the phase is consistent with zero 
as expected in the SM. 

Physics beyond the SM may affect B®B® mixing and 
CP violation in B — > J/ipK® and B — > tttt. Using a model- 
independent analysis [T2j we can introduce a scale param- 
eter, 7*d, for B®B® mixing and an additional phase, 0d, for 
parameterizing a^°- I n the SM we expect rd = 1 and 
0d = 0. With present uncertainties and 0d are consis- 
tent with the SM expectations (see [TO]). 

5 Visualizing the role of theoretical errors 

In addition to the global fits in the p — fj plane, we explore 
the impact of measurements on the theoretical parameters 
and their correlations. We typically scan theory parame- 
ters within ±5Z\ and denote them with ~. Presently, we 
use either exclusive or inclusive V u b , V c b information and 
plot contours for three of the five scanned theoretical pa- 
rameters for different conditions. An example is shown 
in Figure 01 where we have scanned inclusive V u b, inclu- 
sive V cb , Bk, fB d \J BB d and f. For V ub , V cb and B K we 
plot two-dimensional contours on the surface of a cube. In 
each plane five contours are visible. The outermost con- 
tour (solid black) results from requiring a fit probability 
of > 32%. The next contour (also solid black) is obtained 
by restricting all other undisplayed theory parameters to 
their allowed range of ±1Z\. The third solid line results 
by fixing the parameter orthogonal to plane to the al- 
lowed range, while the outer dashed line is found if the 
latter parameter is fixed to its central value. The inter- 
nal dashed black line is obtained by fixing all undisplayed 




Fig. 3. Contours of the theory parameters Bk — V u b — V c b both 
resulting from inclusive reduced semileptonic rates for fit prob- 
abilities P(x 2 ) > 32% after scanning B k Jb^Bb^P[B -> 
X c tv and F{B — > X u iv over ±5Z\* range. 

parameters to their central values. Further details, other 
combination plots and results for exclusive V ub and V cb 
scans are discussed in [TU] . 

6 Conclusion 

The scan method provides a conservative, robust method 
that treats non-gaussian theoretical uncertainties in an 
unbiased way. This reduces conflicts with the SM result- 
ing from unwarranted assumptions concerning the theo- 
retical uncertainties, which is important in searches for 
new physics. The scan methods yields significantly larger 
ranges for the p — fj plane than the Bayesian method. 
Presently, all measurements are consistent with the SM 
expectation due to the large theoretical uncertainties. The 
deviation of a^^o from sin 2(3 measured in charmonium 
Kg (if£) modes is interesting but not yet significant. Model- 
independent parameterizations will become important in 
the future when theory errors are further reduced. 
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